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THREE-DIMENSIONAL STRUCTURE OF SOLAR WIND TURBULENCE 
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ABSTRACT 

We present a measurement of the scale-dependent, three-dimensional structure of the magnetic field fluctu- 
ations in inertial range solar wind turbulence with respect to a local, physically motivated coordinate system. 
The Alfvenic fluctuations are three-dimensionally anisotropic, with the sense of this anisotropy varying from 
large to small scales. At the outer scale, the magnetic field correlations are longest in the local fluctuation 
direction, consistent with Alfven waves. At the proton gyroscale, they are longest along the local mean field 
direction and shortest in the direction perpendicular to the local mean field and the local field fluctuation. The 
compressive fluctuations are highly elongated along the local mean field direction, although axially symmetric 
perpendicular to it. Their large anisotropy may explain why they are not heavily damped in the solar wind. 
Subject headings: magnetic fields — MHD — plasmas — solar wind — turbulence 



1. INTRODUCTION 



The solar wind is a weakly collisional plasma (e.g., 
Kaspe ret al.l |2008|) that is ubiquitously observed to be in 
a turbulent st ate (iTu & Marschll l995t [Goldstein et all 1 1995b 



Horbu rv et alJl2005b [ Bruno & Carbone 2005: iPetrosyan et alJ 
20 lOt iMatthaeus & Vellil 1201 lh . Much progress has been 
made in understanding the nature of this t urbulence since the 
first direct spac ecraft observations (e.g., ISiscoe et al.l 1 19681 : 
IColeman1ll968l) but many aspects remain to be fully under- 
stood. In particular, the three-dimensional (3D) structure 
has been poorly characterized. Here, we use a new single- 
spacecraft technique to measure the 3D structure of turbu- 
lence in the fast solar wind. 

Turbulence is usually modeled as a local cascade of fluctu- 
ations from large to small scales, forming an inertial range. 
In the solar wind, mos t of the energy at large scales is 
in Alfvenic fluctuations ([Belcher & pavis|ll9TiT :TBruno et al. 
ll985UHorbury et alJl995blBale et alJ2005l) . which have mag- 
netic field and velocit y fluctuations p erpendicular to the mag- 
netic field direction dAlfvenl Il942l) . Early isot ropic mag- 
netohydrodynamic (M HD) turbulence theories (Iroshnikov 
119631: Kraichnan|| 1965l) ba sed on Kolmogorov scaling argu- 
ments (Kolmogorov 1941) predict that the energy spectrum 
of weak Alfvenic turbulence is E(k) ~ £~ 3 / 2 , where k is 
the wavenumber of the fluctuations. Although ID velocity 
pow er spectra in the solar wind at 1 AU display this scal- 
ing dMangenev et al J 120011: iPodesta et al J 120071: ISalem et alJ 
120091: IChen et al.ll201 lat iBoldvrev et al.ll201 llCfhe magnetic 
field has a fc~ 5/3 scaling (e.g.. IMatthaeus & Goldstein! 1 19821: 
ISmith et al.l l2006tlChen etaLl l201 latlBoldvrev et alJl201 Hi- 
lt was later re alized ([Montgomery & Tur ned 119811: 
IShebalin et al.l Q983) that the magnetic field direction can 
induce anisotro p y in plasma turbulenc e . It wa s then proposed 
dffigdonl [19841: iGoldreich & Sridhari fl995l) that Alfvenic 
turbulence tends towards a state of critical balance, in which 
the timescale of the Alfvenic fluctuations propagating along 
the magnetic field is equal to the timescale of their nonlinear 
decay. This produces a spectrum perpendicular to the local 
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Figure 1. Normalized B-trace structure function in three orthogonal direc- 
tions. The grey dashed lines indicate the range of values over which spectral 
indices were fitted. The red dotted lines correspond to the eddy shapes in 
Fig.l] 

—9 2/3 

E(k\\) ~ Km and local wavevector scaling k\\ <~ kj_ . Solar 
wind turbulence measurements show evide nce for both 
wavevector anisotropy of the form kj_ > k \\ (ICrooker et alJ 
1 19821: iBieber et alj|1996t iLeamon et all 1 199a Borburv et all 
200a iPodestal 120091: IWicks et all 1201 a IChen et al.l 1201 lbt 
I Wicks et al.ll201 It iHorbury et al.112.01 ll) and a s teeper spectral 
index parallel to the local mag netic fi eld (IHorbury et all 

i et aUboioE 



IPodestal 120091 ILuo & Wul l2010t IWicks ! 



IChen etal]l201 lbt IWicks et al1 l201 lhlHorburv et alJl201 HT 

Critical balance theory was later extended to allow for 
the possibility that Alfvenic turbulence is 3D anisotropic 
(Boldyrev 2006). The two special orthogonal directions are 
the mean magnetic field Bo and the perpendicular magnetic 
field fluctuation <5B^. The theory assumes that the magnetic 
field and velocity fluctuations align to within a scale depen- 
dent angle v b, which makes them 3D anisotropic: / > £ > A, 
where I, £ and A are their correlation lengths in the mean 
field direction Bo, in the 5B± direction and perpendicular to 
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Figure 2. Normalized B-trace structure function at 1.5 X 10 Hz as a func- 
tion of 6b and 8$b ■ 



Figure 3. B-trace spectral index between normalized structure function val- 
ues of 0.08 and 0.3 as a function of 9b and Bgg . 



both, respectively. The local spectra implied by the theory 



in these three directions are E(ki) 



E(ks) 
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The ^-direction scaling follows from substi- 



E{k x ) 

tuting £ ~ A 3/4 into Sv ~ A 1/4 from lBoldvrevl d2006l) . giving 
5v ~ £ 1//3 , corresponding to a local -5/3 spectrum. 

Scale dependent alignment has been reported in the solar 
wind at large scales but is difficult to measur e deep in the 
inertial range due to instrumental limitations (Podesta et al. 
2009). A recent multi-spacecraft measurement of the turbu- 
lent energy distribution in the near-Earth solar wind suggested 
that there was anisotropy with respect to global directions of 
the system, suc h as the global mea n field, solar wind flow or 
the bow shock (Narita et al. 2010b c). As far as we are aware, 
there has not yet been a measurement of the 3D structure of 
solar wind turbulence in a local, scale-dependent coordinate 
system (7,£, A). 

Although inertial range solar wind turbulence is predom- 
inantly Alfvenic, there is also a non-negligible spectrum of 
compressive fluctuat ions SBh and 6n, where n is the num- 
ber density (e.g.. iMarsch & Tul IT990fc ITu & Marschl 11994 
Bavassano et al J 1120041 iHnat et all 120051: iKellogg & Horburvl 
2005t llssautier et al.l 120101: IChen et afTl201 laU2012l). The 
nature of these fluctuations is debated (Matthaeus et al. 



199TtlLithwick & Go ldreich 2001; Kellogg & Horburv i2005T 
Schekochihin et all 120091: iHowes et all 120121: iKleinetall 
20121) . in particular, the reason why they are not heavily 



damped. Their structure has been less comprehensively char- 
acterized than the Alfvenic turbulence, although measure- 
ments in the mag netosheath show that there is some deg ree 
of 2D anisotropy dAlexandrova et al.ll2008h iHe et alj|2011h . 

In this paper, we present measurements of the scale- 
dependent 3D structure of the Alfvenic and compressive mag- 
netic field fluctuations with respect to a new local coordinate 
system and discuss the implications for our understanding of 
plasma turbulence. 

2. METHOD 

In t he analysis, fast sol ar wind data from the Ulysses space- 
craft {yVenzel et al. 1992) during a polar pass between 1 .4 and 
2.6 AU in days 10 0-299 of 1995 wa s used. The magnetic field 
data from VHM dBalogh et al Jl 1992b was at 1 sec reso lution 
and the velocity data from SWOOPS dBame et alJI 19921) was 



at 4 min resolution. The average solar wind speed was ~ 780 
km s" 1 and the outer scale cross-helicity was moderately high, 
a c ~ 0.6 (other pla sma parameters for this stream are given in 



(otner pi; 

iWicks et al.1 d2010l) ). The data was split into 10 equal length 
intervals for the analysis. 

For each 20 day interval, 21 logarithmically spaced 
spacecraft-frame frequencies at which to measure the power 
levels, over the range 5 x 10~ 6 Hz < / sc < 5 x 10 _1 Hz, were 
chosen. For each of these frequencies, the pairs of magnetic 
field measurements, Bj and B2, with the time lag l// sc were 
selected. For each pair, the contribution to the second or- 

der B-trace structure function ^ (Bij-B2j) , where i is the 

i 

component of the magnetic field, and the contribution to the 

second order |B| structure function (|Bi | — | B2 1 ) were calcu- 
lated. Since most of the energy is in the perpendicular fluc- 
tuations (e.g.. iBelcher & Davis! 11971b . the B-trace spectrum 
is a good proxy for the Alfvenic <5Bj_ spectrum and since 
|B| = |B + (5B| w v /|Bo| 2 +2B -5Bw|B |+5B||, the |B| spec- 
trum is a good proxy for the compressive <5Z?|| spectrum in the 
inertial range, where |<5B| < \Bq\. 

It has been shown that using a local scale-dependent 
2D coordinate system is important for testing theoretical 
predictions of spec tral an isotropy dCho & Vishniacl 120001: 
Maron & Goldreichl l200lt iHorbury et alJ 12008b IChen et alJ 
2010a, 1201 lbb . Here, this is extended by defining a lo 



cal scale-dependent 3D coordinate system. For each pair of 
points, the local mean field Bi oca i = (Bj +B2) /2 and the local 
perpendicular fluctuation direction Bi oca i x [(Bi - B2) x Bi oca i] 
were calculated. The angle between Bi oca i and the mean so- 
lar wind velocity (which is the sampling direction), 0b, and 
the angle between the local perpendicular fluctuation and the 
component of the solar wind velocity perpendicular to Bi oca i, 
05B ± , were then found. 

An orthogonal spherical polar coordinate system was de- 
fined, in which / sc is the radial coordinate, 8 B is the polar 
angle and 6sb± is the azimuthal angle. It is in this local 
coordinate system that the 3D anisotropy of the turbulence 
was measured. The structure function contributions for each 
/ sc were binned in 10° bins of 9g and 6sb ± and the mean 
value in each bin was calculated. While the structure func- 
tions conditioned to the local field direction in this way may 
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Figure 4. Surfaces of constant B-trace power (statistical Alfvenic eddy shapes) from large (left) to small (right) scales, in which color represents distance from 
the origin. The normalized power levels are 0.7, 0.17 and 0.02 as marked with red dotted lines on Fig.[T] The typical proton gyroradius is 360 km. 



not be purely second order (Matth aeus et al.ll2012l) . they are 
thought to be the rele vant q uantities in critical balance the- 
ories (ICho & Vishniadl2000t iHorbury ei"aTll20()! loien et al l 
1201 lbl) . Any angles greater than 90° were reflected below 
90° to improve accuracy for scaling measurements. Reflec- 
tion in 9sb ± was found to be a good approximation; while 
there were few points to check the validity of refl ection in 8b, 
the assumption seems reasonable (Podesta 2009). 

Taylor's hypothesis (iTaylorll 19381) can be assumed for this 
analysis: since the speed of the solar wind moving past the 
spacec raft is more than 1 times the Alfven speed in this in- 
terval (IWicks et al.ll2010l) . temporal variations measured by 
the spacecraft, l// sc , correspond to spatial variations in the 
plasma v sw // sc , where v sw is the solar w ind speed. This ha s 
been shown to be a good approximation dNarita et al.ll2010ah . 

3. RESULTS 
3.1. Alfvenic Fluctuations 

Fig. Q] shows the B-trace structure function ("power") as 
a function of spaceraft-frame frequency for three angle bins 
corresponding to the Bi oca ] direction (red circles), the SB± di- 
rection (green diamonds) and the direction perpendicular to 
both (blue squares). Each value is the mean calculated from 
the 10 intervals and the error bars are 2a, where a is the stan- 
dard error of the mean. Before averaging, the structure func- 
tions of each interval were normalized to the square of the 
mean field strength over the interval (|B|) 2 to account for the 
varying power levels due to the spacecraft orbit. The typical 
proton gyroscale corresponds to a spacecraft-frame frequency 
» 0.3 Hz. 

The perpendicular (blue) curve is characteristic of fast 
solar wind: shallow in the low frequency 1 / / sc range 
(Matthaeus & Goldstein 1986) and steeper in the higher fre- 
quency inertial range. The parallel (red) curve also matches 
previous parallel spectrum measurements, following the per- 
pendicular curve at low frequ encies, then becom ing steeper 
than it in the inertial range ( IWicks etalJl20Toh . The <5B ± 
(green) curve has not previously been measured and describes 
how the 3D anisotropy evolves in the turbulent cascade. At 
large scales it has a smaller value than the other structure func- 
tions, which is consistent with this range consisting of Alfven 
waves dBelcher & Davis|[l97ll) . since they have wavevectors 
in the plane perpendicular to 5K±, It also remains smaller 
than the perpendicular structure function throughout the cas- 



cade but becomes larger than the parallel one at « 3 x 10 
Hz. 

For each 20 day interval, a power law was fitted to the nor- 
malized structure functions between values of 0.08 and 0.3 
(marked as grey dashed lines) in each angle bin. A fixed 
power range, rather than a fixed / sc range, was used so that 
the scaling was m easured for the same set of fluctuations 
dChen et al.l 1201 Obi) . For each angle bin, the fit to the struc- 
ture function was evaluated at 1.5 x 10~ 2 Hz to give the 3D 
power anisotropy and the mean of the 10 intervals is shown 
in Fig. |2] The typical standard error of the log of the mean is 
between 0.05 and 0.07. It can be seen that the power increases 
with both 8b and 9sb ± , indicating 3D anisotropy, and seems 
to peak near 8 B = 60°, 6 SB ± =90°. 

Each fitted power la w index was converted to a spectral 
index by subtracting 1 (Monin & Yaglom| [l975l) and the 3D 
spectral index anisotropy is shown in Fig. [3] The typical stan- 
dard error of the mean is 0.01 or 0.02, although the actual 
uncertainty may be larger due to systematic effects, such as 
the finite frequency respons e of the structure fun ctions. The 
steepening towards small 8b (Horb urv et aU2 008) can be seen 
but there appears to be little variation with 6$b ± at large 8b- 

To visualize how the 3D anisotropy varies through the tur- 
bulent cascade, surfaces of constant power were calculated. 
At a selected structure function value, the corresponding fre- 
quency in each angle bin was found through linear interpola- 
tion and the scales corresponding to these frequencies were 
calculated using Taylor's hypothesis. The scales, together 
with the angles 9b and 9$ B± , were converted into Cartesian co- 
ordinates (/, £, A) and the surfaces of constant power (at struc- 
ture function values marked by red dotted lines in Fig. [TJ are 
shown in Fig. [4] They have been reflected into the other seven 
octants under the assumption of reflectional symmetry (see 
earlier). These statistical surfaces can loosely be considered 
as average eddy shapes (although they are not eddies in the 
dynamical sense). It can be seen that they change from being 
extended in the <5Bj^ direction in the large scale Alfven wave 
range (£ > l,X) to being 3D anisotropic close to the proton 
gyroscale (/ > £ > A). 

3.2. Compressive Fluctuations 

The results of a similar analysis for |B| are shown in Fig. [5] 
(where the structure functions have been normalized in the 
same way as for the Alfvenic fluctuations in Fig. Q]) and 
Fig. |6l which is the surface of constant normalized power of 
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Figure 5. Normalized |B| structure function in three orthogonal directions. 
The grey dashed lines show the range of values over which spectral indices 
were fitted. The red dotted line corresponds to the eddy shape in Fig. [6] 

1.2 x 1(T 3 (marked as a red dotted line in Fig. [5j. It can be 
seen that the structure of the compressive fluctuations is dif- 
ferent to the Alfvenic fluctuations: there is no anisotropy in 
the plane perpendicular to the mean field, meaning that the 
compressive fluctuations do not depend on the polarization 
of the Alfvenic fluctuations. Also, they are more elongated 
along the mean field direction than the Alfvenic fluctuations: 
for a given perpendicular scale A, the ratio l/X is at least 2 
or 3 times larger. Due to limited angular resolution this is a 
lower limit; by extrapolating the shape in Fig. [6] one could 
imagine that they are even more extended than can currently 
be measured. 

The spectral indices of |B| for normalized powers between 
3 x 10~ 3 and 1.1 x 10" 2 are between -1.58 and -1.42 in all 
angle bins, with a typical standard error of the mean of 0.02. 
This is different to the slow solar wind, where spec tral in- 
dicies close to -5/3 are observed (IChen et al.ll201 lal) . This 
differ ence has also been n oticed in the electron density spec- 
trum (llssautier et al.ll20Toh . although the reason is not well 
understood. If the compressive fluctuations are indeed very 
anisotropic, then we would not expect to measure the true par- 
allel spectral index with the current angular resolution, which 
may explain the presence of anisotropic structures yet no sig- 
nificant anisotropic scaling. 

4. DISCUSSION 

We have shown that the Alfvenic turbulence is locally 
anisotropic in the plane perpendicular to the mean field. Since 
the direction of the anisotropy is associated with <5B^, the 
question naturally arises to what extent this anisot ropy is a re- 
flectio n of the solenoidality of the magnetic field ( Turn er et al.1 
1201 lb . While the magnetic field has zero divergence at each 
point, this does not imply that the correlation length along 
SB± is infinite, because at any given scale we are considering 
finite field increments, not derivatives. 

The effect of solenoidality can be derived from knowledge 
of the probability density function (PDF) of 5R±. Since we 
are considering the conditional structure function at each scale 
{SB 2 ± \9sb ± ), the joint PDF p(5Bj_,9sb ± ), or, equivalently, the 
PDF of the vector SB± is required. As a simple illustration, 
consider the case where this PDF is Gaussian and, therefore, 
fully determined by the second-order longitudinal correlation 



Figure 6. Surface of constant |B| power (statistical compressive eddy shape) 
at small scales. The normalized power level is 1.2 X 10 -3 as marked with a 
red dotted line on Fig. [5] 

function Czx(r) = (\5B± ■ */r\ 2 ), where r is the point sepa- 
ration in the perpendicular plane. The conditional structure 
function becomes 



(5B\{f)\6 SBx ) 



2C LL (r)C T Ar) 



C LL (r) sin 2 6 SB± +C TT (r) cos 2 6 SB± 



(1) 



where the transverse corr elation function is Ctt(/) = 
[rCii(r)]' from solenoidality (Batc helorll953l) . If, in the iner- 
tial range, Cll(/) oc r a , then 



(SB\(r 



1 + a cos 2 ( 



(2) 



Therefore, the ratio of the correlation scales along and across 
the fluctuation direction for a given structure function value 
is scale-independent and equal to £/A = (1 + a) 1 '", which, 
for the measured value of a ~ 2/3, gives an anisotropy 
in the perpendicular plane of » 2.15. Since the measured 
anisotropy is larger than this (varying between 3.2 and 3.8), 
non-Gaussianity is required to explain the observations. Any 
scale-dep e ndent alignment, e.g., the dynamical alignment of 
Boldyrev (2006), is likely to require non-Gaussianity and 
therefore be closely related to the intermittency of the turbu- 
lence. 

Our results show some important differences to an ear- 
lier study, which suggested that the small scale fluctuations 
are longest-correlated in one of the perpendicular directions 
and that the spectra l index is different in all three directions 
dNarita et a l. 2010b). Possible reasons for this include the dif- 
ferent coordinate system used (g lobal rather t han loc al), the 
presence of foreshock effects in iNarita et al . (2010b) or the 
different assumptions in the measurement technique. 

The fact that the compressive fluctuations are very elon- 
gated is consistent with the prediction, based on gyroki- 
netic theory, that they are passive to the Alfvenic fluctua- 
tions, but have no parallel cascade along the exact magnetic 
field lines (ISchekochihin et alll2009l) . This may explain why 
there is a compressive cascade in the solar wind: the com- 
pressive fluctuations are expected to be damped a t a rate 
proportional to their parallel wavenumber 7 ~ h \ (Barnes 
[19661: iSchekochMnetail 120091: Klein et all I20T2I) but if k\\ 
is very small then they are not heavily damped and can cas- 
cade nonlinearly. An alternative explanation is that the less 
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anisotropic compressive fluctuations are generated but are 
quickly damped, leaving the highly elongated structures to be 
observed. 
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